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ABSTRACT 


An experimental investigation was conducted to measure soot properties and 
concentration, and infrared signature of the plume of a kerosene/gaseous oxygen rocket 
motor. Multiple-wavelength measurements were used to obtain soot parameters. It was 
found that the multiple-wavelength technique appeared to work well for obtaining soot 
properties and concentration. The results implied that the soot particulate was very porous 
for the most fuel-rich conditions and became more compact as the O/F ratio increased. 
Soot mean diameter (D32) was less than 0.1 um when the O/F ratio was less than 0.75, 
but increased to 0.2—>0.36 um for higher O/F ratios. The soot particle diameters obtained 
from the Scanning Electronic Microscope (SEM) photomicrographs of collected soot 
samples were in qualitative agreement with the diameters measured optically. Soot 
concentration in the plume was found to rapidly increase and then decrease as the O/F 
ratio was increased from 0.3 to 1.1, with a peak near O/F = 0.76. Infrared measurements 
in the 3.5-5 um band indicated a probable shift in the dominant radiation source from gas 


to particulate as the soot concentration increased. 
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I. INTRODUCTION 


The use of liquid propellants for launch-vehicle propulsion systems has generally 
been the preferable choice due to their high specific impulse and the ease for providing 
throttling of the thrust level. The liquid propellant has to be sprayed into the combustion 
chamber, in which combustion efficiency is very dependent on intimate mixing of fuel 
and oxidant. Besides that, it is known that mixtures richer than stoichiometric will give 
best results in terms of specific impulse, because of the more rapid reduction in product 
molecular weight than the combustion temperature. Studies on spray combustion and its 
associated processes have helped designers provide for more efficient and stable 
combustors and to better understand the formation of several so-called pollutants such 
as soot, unburned fuel, nitrogen oxides (NOx) and carbon monoxide (CQ). [Ref. 1] Both 
fuel rich conditions and/or imperfect mixing lead to undesirable residue generation, 
among which soot is often the element of concern since it is often associated with most of 
the infrared emissions from a nie or gas turbine plume. Infrared emiting plumes are 
very undesirable, especially for weapon purposes (missiles and even fighter turbines), 
since they make it easier to detect the existence and/or location of the threat, thus 
facilitating the counter-attack. 

The development of reliable computer codes, which can help propulsion system 
designers predict the amount of soot to be generated, is of high importance in looking for 
the least sooty conditions, as well as for analyzing the use of chemical additives in the 


fuel for the same purpose. These codes generally involve the use of some “reduced 


kinetics” scheme together with more standard flow codes for the combustor and exhaust 
nozzle. To validate these codes, it 1s necessary to provide accurate experimental data on 
soot properties and concetrations over a wide range of combustor and nozzle operating 
conditions. 

The size and concentration of soot particles from exhaust plumes have been found 
to significantly influence their heat signatures. Among the several methods which have 
been tried to accurately measure soot concentration in exhaust plumes, such as isokinetic 
sampling, single particle counters and light scattering, “light extinction” is likely the 
easiest and most reliable within a reasonable range of particle concentrations. 

Cashdollar, et al [Ref. 2] gave a practical procedure for using light extinction for 
this purpose, even though the use of their three-wavelength light extinction technique 
could not provide a unique solution for the soot complex index of refraction, mean 
particle size and standard deviation of an assumed log-normal distribution. Four 
independent parameters require the use of a minimum of five wavelengths. Swenson 
[Ref. 3] used a five-wavelength light extinction technique with reasonable success for 
kerosene and gaseous oxygen. To some extent, the present work is the completion of 
Swenson’s effort in establishing a full procedure to determine soot concentration and 
S1Ze. 

This thesis effort has been aimed at developing and utilizing a reliable and 
relatively easy procedure for soot particle size and concentration determination in the 
exhaust plume of a small kerosene-gaseous oxygen rocket motor. The data were to be 
used by other investigators for validation of numerical prediction codes. 
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Il. THEORY AND BACKGROUND 


A. SOOT FORMATION 


The combustion of liquid sprays occurs in liquid propellant rocket motors, gas 
turbines, diesel engines, industrial furnaces, etc. In all these situations, the formation of 
pollutant residues is a reality and soot is among the most significant of them. 


Soot is a carbonaceous solid produced in pyrolisis and combustion 
systems when conditions are such as to allow gas-phase condensation 
reactions of the fuel and its decomposition products to compete with 
further decompositon and oxidation. .... The formation and emission of 
soot are subjects of some concern for a variety of reasons. Soot particles 
are strongly absorbing and, within a combustor, can enhance 
significantly radiative heat transfer. .... The internal structure of soot 
particles has been examined by high-resolution phase-contrast electron 
microscopy. Near the edge of the particle, bent carbon layers follow 
the shape of the particle surface. Inside the particle, lattice structures seem 
to be located more or less regularly around certain centers between which 
the structure is less ordered. ...The density of the particles may be less 
than 2 g/cm3 due to large interplanar spacings. [Ref. 4] 


The process of soot formation involves several phases, such as “particle inception’ 
(condensed phase material from the fuel-molecules via their breakdown), leading to 
‘surface growth’ and ‘coagulation’. While surface growth doesn’t affect the number of 
particles, coagulation influences the establishment of a particle size distribution. 

In Kuo [Ref. 1] we read that “it has been found from various investigations that 
conditions most favorable for soot formation occur when fuel-rich zones have strong 


temperature gradients.” He also states that “the basic method of reducing soot formation 


in spray combustion is not only to reduce temperature gradients in fuel-rich zones but 


also to reduce the size of zones where strong temperature gradients and fuel-rich 
concentration can arise.” One can see that the flow properties should be very well 
determined and controlled to achieve minimum soot formation. 

In rocket motors, liquid propellant atomization is often accomplished by the 
impingement of liquid streams which exit from small holes in the injector face. The size 
of the holes, the injection velocity and the liquid surface tension all influence the 
resulting droplet sizes. 

“The rate at which the droplet evaporates and burns is generally considered to be 
determined by the rate of heat transfer from the flame front to the fuel droplet surface.” 
[Ref. 5] Just as in the case of gaseous fuel diffusion flames, the burning of condensed 
phases often is controlled by rates of heat transfer and diffusion. 

The luminosity of many flames comes from the presence of carbon particles as 
they radiate at the high combustion gas temperatures. Condensed phase particles present a 
very high emissivity when compared to gaseous combustion products. 

It is commonly known that the first step for high combustion efficiency and low 
soot production is getting the reagents to be well mixed, which 1s better attained with 
turbulent mixing. Assuming that reasonably good mixing is achieved, soot production 
may still arise due to other variables, such as the difference between the chamber wall 
temperature and the chamber core temperature together with locally fuel-rich mixtures 
which may lead to unburned fuel thermal breakdown. 

“Particle size and shape depend on the conditions of the flame. They may be 
almost spherical or in the form of filaments.”[Ref. 6] The radii of these particies 
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range from 5 to 300 nm. Such significant variations in particle shape and size led to some 
of the difficulties faced by investigations of particle sizing techniques. 

The words “carbon” and “soot” have been used by investigators to describe a 
wide variety of solid materials, many of which may contain other elements than pure 
carbon. Since the properties of solids change appreciably with the conditions of 
formation, a wide variety of soot properties has been reported. Soot concentration 
measurements based upon light transmission depend on soot density, a property not easy 
to determine. 

B. MULTIPLE WAVELENGTH LIGHT EXTINCTION TECHNIQUE 
1. Particle Size Distribution 
Particle size 1s usually represented by the diameter (in microns) of a spherical 
particle of equivalent mass. Mie theory is based upon spherical particles and it 
works reasonably well if the particles are not much different from the spherical 

shape. For a collection of particles, size information is often given in the form of a 

frequency distribution plotted against particle diameter. In addition, one may 

define a “probability density function - P(D) ” which represents statistically the 
distribution. 

To characterize particle size distributions various diameter definitions may be 

adopted. They are referred to as mean, average, or equivalent diameters. If a 


“number distribution” N(D) is defined as N(D) = N * P(D), where N is the 


total number of particles in the distribution, a general mean diameter for a 


continuous distribution has the form: 


[N(D) D’dD 
2D =. 
JN (D) D*dD 





- (eq. 2.1) 
or for a discrete distribution it would be: 
_2N(P) D'AD 
Pq ya N ( D) IBS AD 
=(Eqe22) 


where N(D) is the number per class width AD [Ref. 7]. Among the possible ‘pq’ 
diameter definitions, the most widely used is the so-called Sauter or Volume-to- 
Surface Mean Diameter (D32). 

From Powell, et al [Ref. 8] one can realize that for a D32 range up to 2.0 
microns “the dependence of the intensity ratio (transmittance) upon D32 is nearly 
independent of the shape of the size distribution function... .” Often it is assumed 
that soot particles have a log-normal distribution. This distribution can be 
characterized by two parameters: the geometric mean diameter (Dg), which 
establishes the central tendency, and the geometric standard deviation (og), which 
determines the spread about the center. 


The log-normal distribution may be presented in the form [Ref. 7 }: 


In(D/ Dg) | 


— — oo =! |e 
one /2 In(og)D _ V2 In(oz) 


ul (4 eo) 


2. Bouguer’s Law 

The best way for obtaining particle size data from in a plume is certainly 
through some nonintrusive method, which in general leads to optical procedures. 
Various researchers [Ref. 9] have made use of the attenuation of monochromatic 
radiation to obtain the properties of particulates. These measurements have been 
usually based upon Bouguer’s Law, also known as the Beer-Lambert Law, which 
states that transmission of light through a cloud of uniform (monodisperse) 


particles is given by: [Ref. 2] 


Pom gi Otned 2. Gees 
- (eq. 2.4) 

where: A - cross-sectional area of the particle 

Cv — - volume concentration of particles 

D - particle diameter 

L - path length containing the particles 

n - number concentration of particles 

Q - dimensionless extinction coefficient 


It should be noted that Bouguer’s Law is defined for monodispersed 
particulate which almost never is the real situation. The extinction coefficient (Q) 
includes loss of light due to both absorption and scattering. 


Dobbins [Ref. 9] revised Bouguer’s Law to apply to a polydisperse system: 


Ow C pel 


i => oe 
2D 32 


sede) 


where Qbar is the average extinction coefficient found for a continuous 


distribution by: 


_JO(D)N(D)D ‘ab 
[N(D)D ‘dD 


bar 


- (eq. 2.6) 
or for a discrete distribution by: 
g. -2Q(P)N(D)D ad 
" ~ N(D)D AD 
- (eq. 2.7) 


Even when the log-normal particle distribution is assumed, the extinction 
coefficient (Q) needs to be calculated for each particle diameter. Mie theory is 
often applied [Ref. 2] to light scattering by a single spherical particle in order to 
obtain the extinction coefficient (Q). Q is a function of wavelength (A) , particle 
diameter (D) and its complex index of refraction (m= a — 1.b). The imaginary part 
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(b) is the absorption index. Cashdollar [Ref. 2] also adds that “even though the Q 
values are calculated on the basis of single scattering, it is valid to use these 
values in Bouguer’s Law for the multiple-scattering case at high concentration of 
particles. The law breaks down only when particles are spaced closer than a few 
particle diameters and interact electromagnetically, thus changing the Q values.” 
3. Mie Scattering Code Application 

Cashdollar, et al assumed the particle index of refraction to be one of the 
values previously reported in the literature. Swenson [Ref. 3] made use of their 
method but let the complex index of refraction vary over all possible values. From 
Cashdollar’s three-wavelength transmission experiments, Swenson had to increase 
to five wavelengths which could provide the four independent ratios required for 
obtaining values for the four variables: index of refraction (m= a — i.b), geometric 
standard deviation (og) and Sauter Diameter (D32). Cashdollar’s Mie code was 
adapted by Swenson to handle five wavelengths. This Mie scattering code 
(Appendix A) was used to generate the Q’s and Qbar’s for the set of five 
wavelengths. The ee were 0.3650um, 0.4047um, 0.4358um, 0.5461 um 
and 0.5770 tm, which represent well defined light emission power peaks from a 
white-light mercury lamp. 

Forming the In-transmittance ratio for two wavelengths using equation 2.5, 


results in: 


In(7U,, ” Orar(AisD32,!,0-¢) 


In Aj , Oj» (A ;,Ds2.tsog) 


(experiment) (theory) 


=(CG.2.6) 


where the transmittances (7) are found experimentally for each wavelength 


through the ratio of light intensities ( J) as follows: 


‘pie I (w/ plume y— L (no ..light ) 


I (no... plume ) Ji (no ..light ) 


=H (op P22) 


As Swenson [Ref. 3] observed, “ ...For five wavelengths, there are ten In- 
transmittances ratios, but only four are independent”... “ required to solve for the 
four unknowns D32 , og, and m= a — i. b .” The “best values” for the four 
variables are found by obtaining the best fit for Equation 2.8 . Therefore, the Mie 
code sweeps through reasonable ranges on all the four variables until it gets the 
best combination, which 1s the one that best falls over a 45-degree line through the 


origin when plotting for all cases: 


aa ie \ys Qbar , 
FS i) Qbar , 
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The actual check used by the Mie code for the best fit is the minimization of 
the sum of the squared deviations (SSD) between the natural log-transmittance 
ratios and the average extinction coefficient ratios. 

Once the best fit was attained for a set of transmittance values all four 
variables were saved by the code, and it was run again for that specific condition 
to obtain the average extinction coefficient. 

For best understanding of the whole method we now assume a fictitious 
experiment for which Table 2.1 expresses the experimental results. Figure 2.1 


shows how the evaluation of the best fit correlated the data. 


Wavelength | 0.3650 | 0.4047 | 0.4536 | 0.5461 | 0.5770 
(micron) 

Transmittance | 0.242 | 0.383 | 0.460 | 0.651 | 0.691 
(0 to 1) | . 


| 
Table 2.1 - Transmittance values from a fictious experiment 
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Figure 2.1 - Mie Code Best Fit Results 
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Once the best fit gave the index of refraction, D32 and og for the case run (e.g. 
Figure 2.1), a second run for this specific case gave the average extinction coefficient as a 
function of the Sauter Mean Diameter (D32), shown at Figures 2.2 and 2.3. 

Picking one of the wavelengths, 0.5770 micron for example, a Qbar of 0.05333 is 
found for the best fit D32 value of 0.076 microns. These value of Qbar and D32, as well as 
the respective wavelength transmittance value (T=0.691) were then used in the Beer- 


Lambert formula (eq. 2.5) to get the desired particle concentration. 
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Figure 2.2 — Qbar vs. D32 for 5 wavelengths (overall) 
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Figure 2.3 — Qbar vs. D32 for 5 wavelengths (detail) 


4. Particle Concentration Determination Method 

After determining the light beam path length through the exit plume, Equation 2.5 
can finally be solved for particle volume concentration (Cv). However, the plume 
code that was used to predict soot levels utilized a particle mass concentration 


defined as: 


= Mass oor 


= Cy pP soot 


mass 
Bas pP gas 


The soot density (psoot) was taken as 1.5 g/cm3 [Ref. 3]. The gas density was 


Cm - (eq. 2.10) 


obtained from the perfect gas law, assuming that the plume pressure (Pe) was 1 
atm. The plume exhaust gas temperature (Te) and molecular weight (Me) were 


obtained at the nozzle exit from a Thermochemical Equilibrium (TEP) code 
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[Ref. 10] using the actual nozzle expansion ratio and measured chamber pressure. 
Te was corrected for the C* efficiency (discussed below). The kerosene fuel was 


asssumed to be similar to RP-1/JP-4 fuel. Thus, Cm was obtained by: 


ex ig 
P.M. 


Cie 


=(eq. 2.11) 
To finally get the ratio between soot (carbon) mass and the total mass coming out 


of the nozzle (soot + gas) we use the definition of Carbon Mass Percentage ( %C(gr) ): 


mass Cm 
0 =| ——__ |x100% = 100% 
Yo C (sr) E } ° [els i 


total 





---- (eq. 2.12) 
5. Technique Requirements and Limitations 
For accurate extinction measurements the detector must view only the 
transmitted light, that is, the detector should never see the forward scattered light 
from relatively large particles. Hodkinson [Ref. 11] suggested that the detector 
half-angle field of view be limited to less than one-tenth the angle of the first 
angular minimum tn the Fraunhofer diffraction pattern: 


wavelength 


< (7.0 
Oy, ( Naar ie 


particle 

- (eq. 2.13) 
Fortunately this was not a hard condition to achieve since the light beam that 
passed through the plume to the detector was collected into a small enough slit 
located relatively far from the plume. 
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Running the Mie code for a best fit using a set of measured transmittance 
values sometimes led to multiple solutions with nearly identical values of SSD. 
Some of these solutions gave unreasonable results for particle index of refraction 
(m). This may be related to the fact that the particle distribution may not be log- 
normal and the particles may not be perfect spheres as Mie theory assumes. When 
a solution resulted in a value for “m’ far from those which have been reported in 


the literature, it was discarded. 
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Ill. EXPERIMENTAL SETUP 


ROCKET MOTOR 


1. Propellant Mass Flow Control 

Propellants were commercial kerosene (liquid fuel) and gaseous oxygen. The 
mass flow control of oxygen was performed through the upstream pressure of a 
sonically choked nozzle. A nozzle throat diameter of 0.0635 inch was chosen to 
cover the whole range needed for the experiments. The sonically choked nozzle 


equation 1s: 





where: m  - Gas Mass Flow Rate 
C, ~ Choked Throat Discharge Coefficient 
P,, ~ Choked Nozzle Upstream Absolute Pressure 


A, ~ Choked Nozzle Throat Cross Sectional Area 


R - Gas Constant for Oxygen 


T., ~ Choked Nozzle Upstream Gas Absolute Temperature 


Ty ) - Function of Specific Heat Ratio for Oxygen 


Wy 


The discharge coefficient ((, ) was experimentally measured as 0.97 using 


an ASME nozzle with a discharge coefficient of 0.98 . Taking ‘y’=1.4, which is 


the room temperature specific heat ratio for oxygen, then f(y)=0.68473, and Eq. 


3.1 can be easily solved for 7 to get the oxygen mass flow rate. 
To make sure the oxygen nozzle is choked, P, _ must be less than 


(0.52828)* P , the sonic choking condition when the specific heat ratio (y) 


upstream 
is 1.4. 

Fuel mass flow rate control was provided using cavitating venturis. These 
devices provide repeatable, stable flow rates that are independent of the 
downstream pressure, enabling gross simplifications of otherwise complex flow 
regulating systems. The fluid medium is accelerated through the venturi throat 
where its vapor pressure is reached, at which point instantaneous ‘flashing’ or 
‘cavitation’ occurs. The vapor is then recondensed in the venturi diffuser with low 
pressure loss and no erosive damage. Flow rate is thereby fixed solely by 
upstream pressure and is unaffected by fluctuations downstream. With a fixed 
feed pressure, cavitating venturis act as ‘hard’ flow regulators, which means flow 
rate is only dependent on upstream pressure and fluid properties. 

To cover the whole range of kerosene mass flow rates needed for the 
experiments, four venturis were used, each with different throat diameter. The 
throat diameters were 0.029”, 0.034”, 0.043”, and 0.052”. Just as the sonic choke, 
the cavitating venturi has a minimum upstream-downstream pressure difference to 
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make sure cavitation happens. This pressure differential (typically a minimum of 
200 psi) was established along with the calibration of each venturi using kerosene 
(Figure 3.1) . These calibrations were performed with the aid of a laboratory 
balance and video recorded timing, so that they gave a very reliable measurement 
of the fuel mass flow rate for each upstream venturi pressure tested. A turbine 
flowmeter was also used in the kerosene line as a check for the cavitating 


venturis. 
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Figure 3.1 - Cavitating Venturi Calibration for Kerosene 
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From upstream pressure measurements taken from fuel and oxygen lines, both 
mass flow rates could be obtained and each (O/F) ratio was accurately determined. 
2. Igniter 

Kerosene 1s not easily ignitable at fuel-rich conditions, even when mixed with 
pure oxygen. This became clear when an ethylene-oxygen torch was used for 
igniting the kerosene-oxygen mixture. To provide a more constant and reliable 
ignition, a hydrogen-oxygen torch was used on all firings. Solenoid valves were 
used to turn the torch gases on and off. Adjusted to a blue flame condition through 
needle valves this torch provided a rapid ignition in almost all the firings. The hot 
gas from the torch flowed through the rocket motor injector head into the 
combustion chamber. Figure 3.2 shows the overall experimental setup used in all 
firings. 

3. Injector Head 

Three similar injector heads were designed and fabricated with one central 
(oxygen) and 8, 6, or 5 peripherical (fuel) orifices. Appendix B shows the rocket 
motor drawings in detail, including the injector head sketch. The purpose of 
different numbers of identical orifices was to have an average fuel droplet size on 
the order of 50um with geometric standard deviation of 2.3 (assumed log-normal 
distribution) for all the fuel mass flow rates required to provide the (O/F) range 


studied. The injector design was driven by keeping its “Penetration Parameter” 
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(Ref. 12] within the optimum range (0.5 to 0.8). With this same intention, the 8- 
hole and the 6-hole injector heads were later turned into 4 and 3-hole injectors by 
welding shut alternating holes. Water was flowed through each injector to insure 


proper impingement of the fuel streams. 
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Figure 3.2 - Firings Experimental Setup 


4. Combustion Chamber Sections 

Three water cooled combustion chamber sections were mounted as shown in 
Figure 3.3 . Detailed drawings of the combustion sections can be found also 
in Appendix B. 

Between the injector head and the first chamber section a thin solid ring was 
used for pressure measurement in the very upstream chamber channel where the 


propellant mixing occurred. An identical ring was used between the second and 
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third chamber sections for a mid-chamber pressure measurement. Between the 
first and second chamber sections a thin solid turbulence ring was used to turn the 
first chamber section into a pre-chamber where the mixing and most of the 
combustion were expected to occur. The actual purpose of the turbulence ring was 


to provide a uniform soot distribution in the exhaust plume. 
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Figure 3.3 - Rocket Motor Overall Geometry 
5. Nozzle 
A water cooled nozzle was designed and fabricated with conical convergent 
and divergent sections. Detailed drawing of this device can be found also in 
Appendix B. All parts of the rocket motor were made of stainless steel 304L. The 
cooling system flowed water independently through the nozzle and each of the 
three chamber sections. The water exits were joined together downstream, where 


water flow rate and temperature were measured for each firing. The measured 
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water flow rate and temperature rise were used to calculate the heat loss from the 

combustor. The heat loss was used to correct the theoretical characteristic velocity 

(C*) obtained from the TEP code. A nozzle discharge coefficient of 0.96 was 

experimentally determined using a choked ASME nozzle with a discharge 

coefficient of 0.98. 
B. WHITE LIGHT SOURCE AND SENSOR SYSTEM 

The white light source used was an ORIEL model 66002 100W Hg lamp with a 
set of 1.5 inch diameter condensing lenses. This Hg lamp provides a set of five very well 
defined power peaks centered at the wavelengths of 365.0 , 404.7 , 435.8 , 546.1 , and 
577.0 nm. The collimated white light beam was reduced to an 0.125 inch diameter by a 
pinhole aperture and passed through the rocket plume. A second pinhole apperture with a 
diameter of 0.0625 inch was used after the beam crossed the plume to prevent the 
measurements from being contaminated by forward scattering. Finally the white light 
beam went through a 25um slit of an ORIEL 77400 1/8 meter spectrograph. The beam 
then spread into its spectrum components and was imaged on an ORIEL Instaspec IV 
CCD detector with 1024x256 pixels on 25um centers. Typically, about 50 center rows 
were binned to minimize any beam steering problems as well as to have a proper average 
power transmission for each wavelength. There was negligible variation in power across 
the 50 diodes. Figure 3.4 shows a normalized transmission plot produced by this system, 
from which the transmittances were calculated. Background correction was possible 


through the system software. 
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The wavelength calibration of the Hg lamp power peaks was performed with 
specially designed Melles Griot optical filters matching the desired wavelengths. In all 
the runs the white light beam crossed the rocket plume at a plane 3 mm away from the 


nozzle exit and almost in all cases through the rocket centerline. 
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Figure 3.4 - Typical Transmission Plot 
c. DATA ACQUISITION SYSTEM 
A Microsoft Visual Basic 5.0 Code was developed to perform the control of the 
whole firing procedure, which included user defined torch time, propellant flow times, 
chamber purge time, and delay between starting torch and starting propellant flow. 
The desktop computer running the Visual Basic Flow Control was coupled with a 
National Instruments BNC-2090 Analog Breakout Box, a PCI-MIO-16E-4 500Ks/s 


Data Acquisition Board and a relay box. This provided the control of all valves as well 
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as the 10Hz digitization of all the analogical device readings, such as pressures, 
temperatures, flow measurements, and laser diode voltage. 
D. INFRARED CAMERA 

An AGEMA 870 Infrared (IR) thermal imaging camera was used to view the 
rocket plume and measure its IR emissions from a point of view perpendicular to the 
rocket axis. The IR camera contained a glass filter with a bandpass of 3.5 to 5 um. 
Images from the IR camera were recorded at 25 Hz and saved into a microcomputer 
harddisk. With an assumed emissivity of the plume, the temperature profile could be 
obtained from the measured irradiance. Conversely, the emissivity of the source can be 
calculated with a known temperature. 
E. HELIUM-NEON LASER 

Besides the white light beam, a Melles Griot 6 mW helium-neon laser beam was 
passed through the rocket plume and directed onto a power meter with a laser line 
interference filter attached. This provided an additional wavelength (632.8 nm) to get 
transmission values if needed for the calculations. On all runs the He-Ne laser beam 
crossed the rocket plume in a plane 3 mm away from the nozzle exit and through the 
rocket plume centerline. 
F, VIDEO CAMERA 

A conventional CCD camera set at a 1/10000" shutter speed was used to record 


the plume image at a standard framing rate of 30 Hz. 
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G. SOOT SAMPLING 

Soot samples were collected with probes connected to a vacuum line. The probes 
utilized 0.24:m paper filters. Two '4-inch diameter tube probes with individual filters were 
positioned 4 feet (fore position) and 20 feet (aft position) downstream from the rocket 
nozzle. Both were oriented toward the nozzle exhaust. The purpose of the soot sampling 
was to obtain qualitative data for comparison with the calculated particle diameters. Soot 
samples were sucessfully obtained only in the last four firings, thus not covering the 
whole O/F range tested. The paper filters were gold plated (after the exposure to soot) 
using a gold plasma and a Scanning Electronic Microscope (SEM) was used to obtain 
photomicrographs. 
H. FUEL DROPLET SIZING 

A Malvern Mastersizer system with 100mm lens was used in an attempt to 
measure the droplet size distribution of the fuel. This system uses a 18mm diameter, 
2mW helium-neon laser beam. For this experiment the injector head was exhausted into 
ambient air and the laser beam was positioned 1.75 inches aft of the injector face, which 


was the gap between the injector face plane and the turbulence ring. 
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IV. EXPERIMENTAL RESULTS 


FIRING RESULTS 
1. Performance Parameters 

More than 70 firings were conducted with the 5, 4, and 3-hole injector heads 
covering an (O/F) range from 0.330 to 1.115, from July through October 1998 . 
Not all the firings were successful. Some presented problems mostly related to 
ignition difficulties when an ethylene-oxygen torch was used and others would 
not ignite when the mixture was too fuel-rich (lower O/F ratios). It was believed 
that the introduction of the ‘turbulence ring’ introduced stagnation conditions in 
the pre-chamber which led twice to combustor failures. These failures partially 
destroyed the rocket motor, requiring new chamber sections and turbulence rings 
to be fabricated. There were also some initial firings when the white light system 
was not yet available ae only the helium-neon laser was used. These runs will 
not be referred to since they were used to test the whole apparatus and firing 
. control hardware and software. Table 4.1 shows the performance parameters from 
the firings. The first digit of the run number was the number of orifices in the 
injector. Theoretical characteristic exhaust velocities were calculated with the 
TEP code including a correction for the heat loss to the water cooling system. 


Experimental C* was calculated as: 


Cuo={ PeAaC./'m, -eqi4 1) 


Ai 


where Pc was taken as pressure downstream of the mixing ring and 7 ; is the total 


mass flow rate. Ignition delays were occasionally encountered but steady state 
conditions were still reached during these rocket firings. Appendix C shows the 
pressure-time traces including upstream oxygen choked nozzle, upstream fuel 
venturi, and chamber pressures and the helium-neon laser transmission diode 


voltages (background corrected). 


O2 Kerosene 
Run | (O/F) | Mass Mass P , : . 
# Ratio Flow Flow amet os +, - Comments 
(Ibm/s) | (Ibmv/s) | (psig) z : 


[5-28 | 0,721 | 0.07043 | 0.0977 | 338 | 4390 | 4256 | 103.1 | Normal Burning 
[5-49 | 0.510 | 0.0774 | 0.1525 | 380 | 3604 | 3963 | 90.9 | Ignition Delay 


Table 4.1 - Performance Parameters 
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Figure 4.1 shows how characteristic velocity efficiency Cem) behaved for 


the O/F ratios tested. The data from Table 4.1 indicate that the 4 and 5-orifice 


injectors had nearly the same performance at the same O/F ratio. For O/F ratios 


lower than 0.8 (’ . decreased significantly. 
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Figure 4.1 - Characteristic Velocity Efficiency vs. O/F Ratios 


2. Performance Parameters 
The transmittances, or the ratios between with-fire and without-fire 
transmissions (with background correction) for all five wavelength peaks 


measured are presented in Table 4.2 for all the above runs. 
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The transmittance results were calculated for the part of every run when the 
transmission values reached a relatively steady value, away from _ the 
influence of the torch ignition. On all the runs the 1/8 inch diameter white light 
beam crossed the rocket plume at its centerline, except for runs 4-10 through 4- 


15, and runs 3-01 through 3-06, when it crossed the plume about 80% and 50% of 





















Run # 





Wavelength (micron) 


Wavelength microns) 


Table 4.2 - Measured Transmittances 








the plume width radially away from the centerline, respectively. The offset was 
employed to increase the transmittance for the plumes with the highest opacities 
in order to be more confident of the applicability of the Beer-Lambert law. For 
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example, compare run 5-29 (O/F=0.744 To.36s = 0.0022) to runs 4-14 (O/F=0.750 
T0.365 = 0.0140) and 4-15 (O/F=0.754 To.36s = 0.0130). The plume width 3mm aft 
from the nozzle exit plane was 16.2mm, while for the 80% and the 50% locations 
the plume widths were calculated as 9.7mm and 14.0mm, respectively. 

On runs 3-04, 3-05, and 3-06 plume afterburning occurred, e.g. exhaust gases 
re-ignited when mixed with ambient air. This significantly contaminated the 
helium-neon laser transmittance, and probably also the white light five-peak 
transmission counts. 

B. SOOT SAMPLING 
Soot sampling was introduced only on the latest runs (runs 3-03 through 3-06). 
Figures 4.2 through 4.12 show the SEM pictures taken of the filter papers. On all SEM 


pictures the small dash on the lower right comer 1s the size scale for 0.5 micron. 
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Figure 4.2 - Run 3-03 - Fore Tube Soot Sample 
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Figure 4.4 - Run 3-04 - Fore Tube Soot Sample 
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Figure 4.5 - Run 3-04 - Fore Tube Soot Sample 


ay 
O@678 15KU @.5u 





Figure 4.6 - Run 3-05 - Fore Tube Soot Sample 
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Figure 4.8 - Run 3-05 - Fore Tube Soot Sample 
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Figure 4.10 - Run 3-06 - Fore Tube Soot Sample 
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Figure 4.11 - Run 3-06 - Fore Tube Soot Sample 
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Figure 4.12 - Run 3-06 - Aft Tube Soot Sample 
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(Ce INFRARED EMISSIONS 

Five infrared images are shown below on Figures 4.13, 4.14, 4.15, 4.16, and 4.17 
for (O/F) ratios 0.338, 0.494, 0.658, 0.895, and 1.115, respectively. For each run the total 
radiance from the plume to about 30 nozzle exit diameters was measured (Table 4.3). In 


addition the nozzle exit temperature was calculated for each run using TEP and corrected 


for ae The plume emissivity was varied in the AGEMA software until the plume exit 


temperature matched the calculated temperature. The resulting axial variation in centerline 
temperature profiles are also shown in Figures 4.13 through 4.17, plotted versus distance in 


nozzle exit diameters downward from the nozzle exit plane. 
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Figure 4.13 - Plume Temperature (K) Distribution for (O/F) = 0.338 
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Figure 4.15 - Plume Temperature (K) Distribution for (O/F) = 0.658 
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Figure 4.16 - Plume Temperature (K) Distribution for (O/F) = 0.895 
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Figure 4.17 - Plume Temperature (K) Distribution for (O/F) = 1.115 
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Table 4.3 - Total Radiance Measurements 












Emissivity (¢) | Total Radiance (W) 


(average) 









D. VIDEO IMAGES 


Figures 4.18, 4.19, and 4.20 show the effect of O/F ratio on the plume opacity. 





Figure 4.18 - Visible Plume Image for O/F=0.510 
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Visible Plume Image for O/F=1.097 
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E. FUEL DROPLET SIZING 

Several attempts were made to obtain the droplet size distributions with the 
Malvern Mastersizer system placed from 1.75 inch to 5 inches away from the injector 
head face. Unfortunately all of them presented transmittance values of about 1%, which 
made the analysis throroughly unreliable. Therefore, no good data could be obtained from 
these experiments. Future efforts will utilize a phase-Doppler particle analyzer (PDPA) to 


measure the spatial variation in the particle size distribution. 
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V. EXPERIMENTAL ANALYSIS 


A. SOOT PROPERTIES 


For each firing the transmittance values were run on the Mie code whose 












Sauter 
Diameter 


Index of Refraction 


results are presented in Table 5.1. The parameter SSD (Sum of Squared Deviations) 
(m=atb.1) 
(a +/- 0.02) (Gg) D32(um) 


Run# | (O/F) 
Ratio 
(b +/- 0.01) (+0.04/-0.16)* | (4/- 0.010) 


S| 0330 | 253-0031 | 1.01 | all | 0.04422 


Table 5.1 - Mie Code Results (* omin= 1.00) 


Log-Norm. 
Sigma 
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defined how well the values for index of refraction, log-normal sigma and 
D32 resulted in the best agreement between InT-ratios and Qbar-ratios. A perfect fit would 
yield SSD = 0. The results of Table 5.1 represent the minimum SSD values that could be 


obtained for each set of transmittance values. In general, values of SSD which 


represented very good fits had orders of magnitude of ro" to] om Some of the SSD 
values were greater and were mostly the ones with greater value for the real part of the 
index of refraction. 

The variation of “a’’ (index of refraction real part) with O/F ratios did not present 
a significant trend. The values were typical of those reported in the literature [Ref. 2,13]. 
Values of “a”>1.95 generally occurred when the data correlation was poor. The 
absorption coefficient (“b’’) was quite low for the most fuel-rich conditions. As the O/F 
ratio was increased the absorption coefficient approached values most often measured for 
soot (0.3 — 0.7) [Ref. 2,13]. The low “effective” values of “b” were also observed by 
Powell & Zinn [Ref. 13] and Santoro, et al [Ref. 14]. They report that the effect was due 
to the particulate consisting of loosely packed soot. Powell & Zinn [Ref. 13] determined 
that for a measured value of “b’=0.075 the fraction of optical mean volume occupied by 
particulate was only 0.184. The current data appear to indicate that very fuel-rich 
conditions result in very porous soot particulate which becomes more dense as the O/F 
ratio is increased. Except for three values, when the O/F ratio was less than 0.75, D32 was 
smaller than 0.14sm. When the O/F ratio was higher than 0.75 D32 was 0.20-0.36um, 


which is more typical of the sizes observed in gas turbine engine plumes. The relatively 
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high uncertainty on og shows that even for the very low values (e.g., 1.01) the actual 
particle distribution may be far from monodisperse. 

The uncertainties presented in Table 5.1 for index of refraction (real and 
imaginary part) and geometric standard deviation (og) were obtained using the Mie code. 
Each variable was individually varied away from the optimal solution until the calculated 
D32 had a significant variation (usually 10%). A similar approach was used for the D32 
uncertainty, using a variation in SSD to three times the best-case value. 

B. SOOT CONCENTRATION VERSUS O/F RATIO 

With the methodology presented in Chapter II, the soot mass concentration was 
calculated. The results are presented in Table 5.2 and plotted in Figure 5.1. 

The achievement of low values of SSD was in most cases relatively easy. 
However, for some transmittance sets, if the ranges for index of refraction (m) that the 
Mie code swept through were not limited to reasonable values, even lower values of SSD 
could be obtained for higher (>3.0) values of the real part of ‘m’. These over-calculated 
values for the real part of ‘m’ were far from those presented. anywhere in the literature for 
soot and were not used. The next best fits were taken, still with very low SSD values. It 
should be noted that %C(gr) was obtained using Psoot=1.5 g/cm3. As pointed out in the 
previous section the particulate density may have been considerably less for the most 
fuel-rich conditions. This would decrease %C(gr) below the values shown in Table 5.3 , 


giving a much steeper rise on %C(gr) in Figure 5.1. 


45 











Run? [OF ratio | %C@) | Run¥ | OF ratio | “Can 
(26 [oar [isis 08 Tos81 15.98 
[143610 
[332 | 0658s 30.97 15 40 
70.750 
P5335 | 0338096) aS 5s — 
sa eae ee ea 
= 
pari i ee i758 | 
Psat 
P5490 
401 
402 













x41 | 0.330 080 [456 
549 | 0510 | 12.61 1.097 

401 | 03575 | 16.22 Lis | 938 | 
P40 [0.566 15.71 } 


Table 5.2 - Soot Mass Concentration vs. O/F Ratios 
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Figure 5.1 - Soot Mass Concentration vs. O/F Mass Flow Ratios 
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Using a simple average of the transmittance for all five wavelengths as a 
parameter, it was observed that when it was lower than 4% (which happened for O/F 
ratios between 0.6 and 0.90, as shown in Figure 5.1) there was considerable scatter in the 
soot concentration values, except for runs 3-01 through 3-03. Outside this O/F ratio range 
soot concentrations were significantly repeatable. Runs 3-04 through 3-06 seeemed to 
have been significantly affected by the plume afterburning. Unfortunately this O/F ratio 
condition could not be repeated in time due to the combustion chamber failure, but the 
transmittance values were very stable. For the most fuel-rich condition tested (O/F ~ 
0.33) it was hard to get good results from the Mie Code, which may have been related to 
a significant portion of unburned fuel mixed with the soot particles in the plume or to the 
probable porosity of the soot particle. These were relatively cold firings and it would be 
likely that condensed fuel droplets could exist in the plume. The particulate analysis 
assumed that only soot particles were present in the plume. 

Figure 4.1 shows C* efficiencies over 100%. This behavior comes possibly from 
the theoretical C* values given by the TEP code, which assumes perfect mixing and 
instantaneous chemical waitin. These assumptions may give lower chamber 
stagnation temperature and C*’s than actually exists because of incomplete combustion. 
It is believed that a significant portion of the fuel could actually be going through the 
rocket motor without any chemical reaction. Thus, for the fuel-rich conditions of this 
investigation, combustion would occur closer to stoichiometric conditions. This would 
increase the chamber pressure and, thus, the experimental C* value above those obtained 
assuming chemical equilibrium of all combustion products. 


AT 


&: UNCERTAINTY ANALYSIS FOR THE SOOT CONCENTRATION 

An uncertainty analysis was conducted on the O/F ratios and on the %C(gr) 
values. O/F ratios obviously depended on the oxygen and fuel mass flow rates. Oxygen 
mass flow rate uncertainties were calculated for a +/- 0.0005” tolerance on the choked 
throat diameter, a +/- 2K accuracy on oxygen temperature, +/- 2 psi on upstream 
pressure, and 0.5% on choked throat discharge coefficient. This resulted in a +/-2.53% 
uncertainty in oxygen mass flow rate values. On the fuel side, accuracies were +/- 2 psi 
on upstream venturi pressure, +/-0.94% for the 0.052’’venturi throat diameter venturi and 
+/-0.65% for all other venturis used. The accuracies resulted in a +/-2.1% uncertainty in 
O/F ratio for runs using the 0.052” venturi and +/-1.8% uncertainty in O/F ratio for runs 
using all other venturis. Figure 5.2 shows the previous plot with these O/F uncertainties 
included. Considering now the uncertainties on %C(gr) one of the runs was picked for 
full analysis. Run 4-02 was selected for its relatively low SSD and for being outside the 
lower-than-4% transmittance O/F range. For the mean transmittance values for each 
wavelength obtained from the Oriel system we originally had a best fit of: 
m=1.70-0.191 , og=1.21, D32=0.1002 micron, SSD=0.0008649, Qbar(577nm)=0.26413 
T(577nm)=0.1500 , %C(gr)= 15.71%. 

The uncertainties in %C(gr) may come from the uncertainties in transmission 
values read by the Oriel system, and from uncertainties of the index of refraction, 
geometric standard deviation and soot mean diameter (expressed in Table 5.1). It 


was found that the transmission readings were very steady and presented a negligible 
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Figure 5.2 - Soot Concentration vs. O/F Ratio with Uncertainties 


uncertainty. The influence on %C(gr) was calculated independently for the uncertainties 
in “a’,”’b’, og and D32. Uncertainty in “a” was +/- 0.02, which gave a +/- 4.80% 
uncertainty in %C(gr). Uncertainty in “b” was +/- 0.01, which gave a +/- 3.55% 
uncertainty in %C(gr). Uncertainty in og was —0.16/+0.04, which gave a +3.12%/-1.23% 
uncertainty in %C(gr). Uncertainty in D32 was +/- 0.010um, which gave a_ +/- 5.6% 
uncertainty in %C(gr). The uncertainty, when considered altogether, may increase above 


the independent uncertainty figures, and the unprobable worst case cennario (all pluses/all 
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minuses) would give a maximum uncertainty on %C(gr) of +/- 16.1%. It would be safe 
to say the overall uncertainty on %C(gr) was probably lower than +/- 5%. As discussed 
above, the variation in soot particulate density would have a one-to-one influence on 
%C(gr). 
D. INFRARED EMISSIONS 

Figures 4.13 through 4.17 showed temperature distributions along the plume 
centerline. Temperatures in the plume ranged approximately from 600K to 2200K. The 
Mach disks resulted in a significant part of the total radiation. It should be noted that the 
Mach disks are highly visible both optically and by the temperature profile, even for 
the most opaque plumes. Table 43 presented the total plume radiance measured to 28 
nozzle exit diameters (De) downward from the nozzle exit plane. Even though the 
centerline temperature distribution seemed to behave as 1f no afterburning was occurring 
for the O/F ratio of 1.115, it was observed in the video image aft of 28 De. The measured 
radiance for O/F=1.115 was expected to be significantly lower than for O/F=0.895 (a 
significant drop in soot concentration), which was not obtained. This was probably due to 
the significant radiation from the downstream afterburning that reached the camera lens. 
If this occurred, the temperature would have been much higher than the value of 1010 K 
used in Table 4.3, and would result in a lower emissivity. The low values of emissivity 
are more typical of those for gases and the higher values to those of particulate, possibly 


indicating the dominant radiation mechanism shifts as the soot concentration increased. 
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E. SOOT SAMPLING 

Figures 4.2 through 4.12 showed the photomicrographs of runs 3-03, 3-04, 3- 
05, and 3-06 which had D32 values of 0.3602um (og=1.01), 0.2001lum (og=1.04), 
0.2815um (og=1.01), and 0.330lum (og=1.01), respectively. These low values of og 
would indicate a nearly monodisperse distribution. However, it has already been shown 
that the uncertainties were high on og, which 1s supported by the photomicrographs. 
Figures 4.2 and 4.3 (run 3-03 — fore tube) showed that spherical particles were a 
significant presence with diameters ranging approximately from 0.1 to 0.4m. It also 
showed some filiform particles and particle agglomerates. Figures 4.4 and 4.5 (run 3-04 — 
fore tube) showed not as many spherical particles, with diameters ranging approximately 
from 0.1 to 0.3um, but a greater amount of irregular but compact particle agglomerates as 
big as 0.75um. Figures 4.6, 4.7, and 4.8 (run 3-05 — fore tube) showed not only a 
significant presence of spherically shaped particles with diameters ranging approximately 
from 0.1 to 0.4m, but also filiform agglomerates of particles within this diameter range. 
Figure 4.9 (run 3-05 — aft tube) showed basically the same particle diameter range except 
that the presence filiform particles and agglomerates was more significant. Figures 4.10 
and 4.11 (run 3-06 — fore tube) showed a significant amount of spherically shaped 
particles with diameters ranging from 0.1 to 0.6um. Not a significant amount of 
agglomerates nor filiform particles were observed. Figure 4.12 (run 3-06 — aft tube) 
showed some particles as big as lum which may not have been soot, but the picture 


quality did not allow a better analysis. 
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It is not known whether the filiform particles and agglomerates existed in the flow 
or were a result of the collection process. The larger size of the spherical particles from 
these runs were in qualitative agreement with the 0.20-40.36 um values of D32 obtained 


optically. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


The multiple-wavelength light extinction technique was generally successful for 
obtaining soot concentration and properties in a plume of a kerosene/gaseous oxygen 
rocket motor over the fuel-rich O/F ratio range where the soot concentration was 
predicted to have a rise-peak-drop behavior. 

The results implied that the soot particulate was very porous for the most fuel-rich 
conditions and became more compact as the O/F ratio increased. Soot mean diameter 
(D32) was less than 0.1 um when the O/F ratio was less than 0.75, but increased to 
0.2->0.36 um for higher O/F ratios. The soot particle diameters obtained from the SEM 
photomicrographs of collected soot samples were in qualitative agreement with the 
diameters measured optically. Soot concentration in the plume was found to rapidly 
increase and then decrease as the O/F ratio was increased from 0.3 to 1.1 , with a peak 
near O/F = 0.76. Infrared ene ments in the 3.5-5 tm band indicated a probable shift 
in the dominant radiation source from gas to particulate as the soot concentration 
increased. 

The turbulence ring combustor design apparently could not withstand the 
conditions imposed by the higher O/F ratios and should be improved or re-designed. The 
Malvern Mastersizer was inadequate for the spray particle sizing. A phase Doppler 
particle analyzer (PDPA), which can typically measure particles down to a minimum 
diameter of 0.5m with a dynamic range (dmax:dmin) of 50:1, could be used with much 


better chances of getting reliable results. 
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The experimental technique developed in this investigation is basically ready to 
go to the next step in the planned study, which would analyze the influence of fuel 
additives on soot generation. For better results it is recommended: to shield the light 
sensors from background radiation such as afterburning using a ducted light beam; to 
possibly increase burning time to get very stable transmission results, away from torch 
gas influence and transients; to increase the white-light lamp power and then use selective 
optical filters to equalize the 5 wavelength power peaks to minimize the influence of 


noise, especially on the ultraviolet range. 
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APPENDIX A. MIE SCATTERING CODE LISTING 


The Mie Scattering Code listed hereafter is a modified version of the miescat4.for 
code used by Swenson [Ref. 3]. The data reduction process (i.e., finding the best 
combination of complex index of refraction, geometric standard deviation and Sauter 
mean diameter for the measured transmittances) was initially modified to be run on a 
microcomputer using Microsoft Fortran Powerstation, but the attempt was not successful 
due to a problem with this software which could not be solved in time. So the version 
presented here was run on the NPS Aeronautics Department’s Silicon Graphics 
Powerstation in a Unix environment. The final version was checked for proper results 
against the original Cashdollar program [Ref. 2]. The code allows iteration starting 
values, step sizes, and the number of loops to be input in the data file for complex index 
of refraction and geometric standard deviation. The difference between the current 
miescat5.for and Swenson’s miescat4.for was that miescat5.for did three loops instead of 
two, sweeping through real and immaginary parts of the index of refraction besides the 
geometric standard deviation. Another new feature was that when all numbers of loops 
were Set to 1 in the data file, the program generated and saved Q vs D, and Qbar vs. D32 
for the specified index of refraction and geometric standard deviation. If any of the 
number of loops was not 1, then the code understood it was a search for the best case, did 
not save Q vs. D nor Qbar vs. D32, but saved all the best cases for each set of index of 
refraction and geometric standard deviation. This made the program run faster and 


generate only those files necessary for each situation. 
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Keka K keke x =€6PROGRAM MiBSGATS.F FOR Se@WAVELENGRHS  * <i ~~ ** 


BY KENNETH CASHDOLLAR, 1976 
REVISED BY THOMAS WELDON, JULY, 1977 
BASED ON PROGRAM MICOEF BY C D LITTON AND Z J FINK 


modified by Sailivinoe 1.Ce SILVA tomobtain, Ehe sest 
index of refraction, log-normal sigma and 
respective D32, for a given set of transmission 
coefficients (from 5 wavelengths) 


QBAR CALCULATES ABSORPTION AND EXTINCTION COEFFICIENTS FOR LIGHT 
TRANSMISSION THROUGH DUST PARTICLES WITH COMPLEX REFRACTIVE INDEX 
AND A LOG-NORMAL SIZE DISTRIBUTION, THEN PLOTS THE RESULTS 


KaeKK RK KKK KK KKK KK Ace INPUT FORMAT Kak KK KKK KKK KK KKK KKK KKK KK KK KKK kk KK 


VARIABLE NAMES AS IN PROGRAM 


ROW 1 Tee, Tl 2g Sa 4: 15 

ROW 2 WAV (1),WAV (2) ,WAV(3),WAV(4),WAV(5) 
ROW 3 NOMR, STEPR 

ROW 4 NUMI, STEPI 

ROW 5 NUMS, STEPS 

ROW 6 RDUST, RMED. 

ROW 7 L265 ,DSPDELDS 

ROW 8 ND, SIGMA, SD32, SDL32 


FORMATS AND VARIABLE DESCRIPTION 
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OM) Ore 


PROGRAM MIESCATS5 


COMPLEX RK, RBUST, RSUST1, D1, DELDI, RA, RINDX (35), Dis; DELIS 
COMPLEX SCOCEPR, ORR, Addn DO hr he Crum S | 

DIMENSION “DPRV1L000) , XJRTLO000 bak (1000) PR (1000) 

REAL X,D, DELD, DS, DELDS, CO, COE, CORP, CEXT ORES, Osea AM, AREL, AMINT 
REAL SD32, SDL32,SIGMAV STEER, STEPL STEPS, RDUSL LR RDUSTILE, WAVE 
REAL WAVL(5),QBARCHK, DQ, DIAM32 

INTEGER M,NPTS1,MINT, K, NOMR, NUMI, NOMS, L 

INTEGER ND, NBEST 

REAL XY (1000) ,WAV(5),SMALLEST 

REAL DP(1000), XK(1000), XJ (1000)7 Dd 1000) 22 e 00) 

REAL RAT51 (1000), RAT52 (1000), RAT53 (1000) , RATS4 (1000) 

REAL RAT41(1000), RAT42 (1000) , RAT43 (1000) 

REAL RAT31 (1000), RAT32 (1000) , RAT21 (1000) 

REAL RATS1A(1000) , RAT52A(1000) , RAT53A(1000) , RATS4A(1000) 
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REAL RAT41A(1000), RAT42A(1000),RAT43A(1000) 

REAL RAT31A(1000), RAT32A(1000) , RAT21A(1000) ,SSD(1000) 
REAL GAMMA, BETA 

COMPEE wes 000) ,B(1000),C (L000) ,RDUST1L 

COMPLE eH, XH (O00) 7 DH (1000) ,E(1000) , DE @O00) 
COMMON/ADQ/DQ (1000, 2) 

COMMON ALLOQ(10000, 4,5) 
COMMON/AVG1/DIAM32 (1000) , WAV(5) 


OPEN(UNIT=1, FILE='mie5Sinp.dat', STATUS='UNKNOWN') 
OPEN (UNIT=2, FILE='QOLnTRatios.dat', STATUS='UNKNOWN') 
OPEN (UNIT=3, FILE='QOvsD.dat', STATUS=' UNKNOWN ' ) 
OPEN (UNIT=4, FILE='QObarD32.dat', STATUS=' UNKNOWN') 
OPEN(UNIT=5, FILE='BestD32.dat', STATUS=' UNKNOWN') 


Ce ee ee ee a ee a a a A ee ee ee ee ee ee ee ee ee ES ee 
e 

Sc ROS ip boweoGr LEX REFRACTIVE INDEX OF BUST PARTICLES,” WLITHSMINUS 

cS RMED TS REAL REFRACTIVE INDEX OF MEDIUM 


peep tl, LOZ, 12,053,164, 15 
HRETE Soe rhe re rs) T4,7T5 
OZ Peeled (ieee a 2h, bso, Z2X,F/7.5,2X, F725, 2X, F7.5) 
READ (1,103)WAV(1),WAV(2),WAV(3),WAV(4),WAV(5) 
WRITE (*,103)WAV(1),WAV(2),WAV(3),WAV(4) ,WAV(5) 
Ios PORMAT (EGe4, 2X,F6.4,2X,F6.4,2xX, F6.4, 2X, F6.4) 
Werte (5, c0liwayv (1), tl, wav(2),t2,wav(3),t3,wav (4) aie4, Wee 5) 7 t5 
READ (1,104)NUMR, STEPR 
Wikre (oO hOs ele ato, F4,5 
write(5,103)wav(l1),wav(2),wav(3),wav(4),wav(5) 
write(5,104)numr,stepr 
WRITE (*,104)NUMR, STEPR 
READ (1,104)NUML,STEPI 
write(5,104)numi,stepi 
WRITE (*,104)NUMI,STEPI 
READ (1,104)NUMS, STEPS 
write(5,104)nums, steps 
WRITE (~, 204) NUMS,STEPS 
104 FORMAT (23, 2X, F4.3) 


READ (1, 32) RDUST, RMED 

write(5,32)rdust,rmed 

WRITE ("7 o72) RDUST, RMED 
a2 FORMAT (2F6.3,2X, F6.3) 


READ 4) b2,ls,Ds, DELDS 
write(5, 34)02,L3,ds,delds 
WR ns 4) b2,83,Ds, DELDS 

34 BORMAT ls, 2%, 15,2X,F0.4, 2X, FC. 4) 


READ (1,12)ND, SIGMA, SD32, SDL32 
Waee( 52 nad, sigma, sd32,sdL32 
WRtte.-) 2) ND, SIGMA, SD3Z,SDL32 

12 BORMAim@legweZ ab 5.3, 2X,F5.3,2X,F5.3) 


RDUST1=RDUST 
RDUST1R=REAL (RDUST1) 
RDUSTLI=IMAG (RDUST1) 


>/ 


SMALLEST=1000. 
OPTSSD=1000. 


et tes crm cer cr cr et et rm cr cr crs crs crc crs cr cr rr cr rr cr rr a a cr cr rs cr cr crs cr cr cr cr me es cs ee we ee ee eee ee 
et rs cr cr cr cr et mc cr cr cr cr cr cr ce ee ee a es ee ee ee eon one es es oe SS ee oe ee Oe eee 


DO 9 INDEX1=1,NUMS 
RDUST1R=REAL (RDUST) 
write (*,655)indexinums 
655 format (/,*ON Sigma loop Index!=",i37e/ ji 
(STARTING OVERALL LOOP ON SIGMA) 
DO 8 INDEX2=1,NUMR 
RDUSTLI=IMAG (RDUST) 
write (*,656) index2, numr 
656 format(' On RauUsER@lceo index2— as, 0/7 mes 
(STARTING OVERALL LOOP ON RDUSTR) 
DO 6 INDEX3=1, NUMI 
write (*,657)index3, numi 
657 focmat ( On RdUSED  leop. badexs—',13,"° /2 70s 
(STARTING OVERALL LOOP ON RDUSTI) 
IMAGINARY PART OF R CHANGED TO POSITIVE TO CONFORM TO EXPANSION 
FORMULA IN PROGRAM 
RDUST1=CMPLX (RDUSTIR, 0.0) +CMPLX (0.0, RDUSTI1T) 
RX=RDUST1/RMED 
R=CONJG (RX) 
RINDX (1) =RX 
RINDX (2) =RX 
RINDX (3) =RX 
RINDX (4) =RX 
RINDX (5) =RX 
INITIAL VALUES FOR ARRAY 
DO} 51 K= 
DO 52 J=1,2 
PO { 1,3, O00 
ALLQ(1,d+2, K) =0.0 
DO 53 I= 2,10000 
ALLO (I,J, K)=100.0 
ALLOW ot2, Kar. 
53 CONTINUE 
52 .CoOnNL TENE 
Si CONTINUE 
DO 54 I=1,1000 
DIAM32 ( 0.0 


erect alli 


li tt CoC OOO Oe @ oc | 


OOO 2 @ ec 


H—_—~_——eE er ess an ss aes ae” 
oOo @ COC oe. © 


Oo Qe © eae 
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(CE GeO) Oa 


DO 471-5 
(STARTING OVERALL LOOP ON WAVELENGTHS) 
WRITE(2,106) Wav (NI) 
WRITE (*,106) Wav (NI) 
106 FORMAT(/,' Entering the Wavelenth loop for Lambda=',F7.4) 


WAVE=WAV (NI) 
WAVL (NI) =WAVE*10000.0 


INITIAL VALUES SET 


MAX VALUE FOR D IS (1000*WAVE)/(1.3*PI*RMED) = 244*WAVE/RMED 
D IS REAL DIAMETER IN MICRONS, DELD IS INCREMENT 

Dl eto COMPLEX PIAMETER, DELDI IS INCREMENT 

X IS REAL SIZE PARAMETER X= PI*DIAMETER*RMED/WAVELENGTH 
ALPH IS COMPLEX SIZE PARAMETER, ALPH=1*X 


D=DS 
DELD=DELDS 
DI=DIS 
DELDI=DELIS 


De Sey, 

Dome f—1 bz 

X=3.14159*D/WAVE*RMED 
ALPH=3.14159*DI/WAVE*RMED 

Xo 1) = Six) 7X 

Mee 2) — oman (0 Ge 2) = COS (XK) /X 
XR} =CSIN (Re) 7 (R*x) 

XJR (2) =CSIN (R*X) / ( (R*X) **2) -CCOS (R*X) / (R*X) 

XH (1) =CEXP (ALPH) /ALPH 

XH (2) =( (-CEXP (ALPH) ) /X)-((0.,1.) *CEXP(ALPH) / (X**2) ) 
Donal) =" COs (x) 7 x— SIN (X) / (X**Z) 

DJR (1) =CCOS (R*X) / (R*X) -CSIN(R*X) / ( (R*X) **2) 
DH(1)=CEXP (ALPH) /X+((0.,1.)*CEXP(ALPH) / (X**2) ) 
OE 1—0 © 

QOSCA=0.0 


Ll = ITERATIONS FOR SPHERICAL FUNCTIONS, L1=1.3*xX 
lS , Se 
Paci tia r@ jas L1—10 
Oe? Te ee 
GO—( FLOAT (1) )/ (2*L+1) 
COE= (FLOAT (L+1))/(2*L+1) 
COEF= (FLOAT (2*L4+1))/X 
COEFR= (FLOAT (2*L+1) )/(R*X) 
XJ (L+2) =COEF*Xd (L+1)-XdJ(L) 
XJR(L+2) =COEFR*XJR(L+1) -—XJR(L) 
XH (L+2) =COEF*XH (L+1) -XH (L) 
Doin iy =CO*XJ (1) -COB*XJ (L+2) 
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Pe ak 


DIR (L+1) =€O*XIJR(L) -COn*xIR (E+2) 
DH GG) =CO-xAH (1) -COn xr Es 2) 
P(L) =X*XJ (L+1) 
PR{(L)=(R*X)*xXgR (Ls!) 
DP (L) =XJ (L+1) +X*DJd (L+1) 
DPR(L)=XJR(L+1)+(R*X) *DJR(L+1) 
BE (L)=X*XH (L+1) 
DE (L) =XH (L+1)+X*DH(L+1) 
A(L)=( (DPRa@L) *P (Djs Gre ie Dee) ) / ( (DPR Ge) B(L))-~(R*PR(L) *DE(L 
I), 

B(L)=( (R*DPR(L) *P(L) )—(PR(L) *DP(L) )) / ( (R*BDERIL) *E (LS PR (LL) *DE(L 
ee 

C(L)=A(L)+B(L) 
GAMMA=CABS (A(L) ) 
BETA=CABS (B(L) ) 
eh i= ( FLOAT (24b41))* (2.7 (Xe 2) * (GAMMA 2 eee 2) 
XY (L)=(FLOAT (2*L+1))*(2./ (%**2) ) * (REAL (CD) ® 
OEXT=QEXT+XY (L) 
OSCA=OSCA+XK (L) 

1 CONTINUE 

QABS=QEXT-QOSCA 
DO (M,1)=D 
DO (M, 2) =OQEXT 


WRITE (2, 22)X%,DQ(M,1) , QABS, DO(M, 2) ,M 
22 FORMAT (9X, 2HA=,F8.3, 4X, 2HD=, F/.3,4X, SHOABS=, F6.3, 4X, SHOEXT=, F6.3, 4 
) 
D AND OEBEAT FOR MONODISPERSE POT INTO PLOTTING ARRAY 
EVERY FOURTH POINT CALCULATED IS PUT INTO ARRAY- ALLO 
XM=M 
AREL=XM/4.0 
MINT=M/4 
AMINT=MINT 


IF (AREL.NE.XMINT) GO TO 60 

Tee Dac sl aZ207"CGO TO 60 

NUMBBR OF POINTS FOR PLOTS 1 AND 3 IS@HPrSi 
NPTS1=MINT 


ALLO (M,1,NI)=D 

ALLO (M, 3,NI) =QEXT 

INCREMENT SIZE PARAMETERS D=DIAMETER 
D=D+DELD 
DI=D1I+DELDI 
M=M+1 

2 CONTINUE 


07-02 OO) Oa aaa 


(end loop on L2) 


DELD=2.0*DELD 
DELDI=2 .0*DELDI 
3S CONTINUE 
(end loop on L3) 
CALL QAVG(M,NI,ND, SIGMA, NPTS2, SD32, SDL32) 
K1=ND+1 
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4 CONTINUE 


c (end loop on NI) 
& cm ce ae ee ee ee ee we me we we we we we we ew ew we we ow ow ow ws ow ow ow we ow ws ow we ow ws ow ws we we = ee we ew we ee ew we we 
S This DO(600) is to print out QvsD , when it is 

the case of running a single value for index of refraction(m) and sigma 
(in this case, m and sigma should be set the wanted value and all the 
variables numr,numi and nums made equal to 1) 
G BO OCOse-17 1-1 
C Weise aGel) “DLO, 1,1) AlebO(r,3,1),ALLO(1, 3,2) ,ALLO(1, 3,3) ,ALLO(I 
e pei CO) 5 GINO M Giles Fie a 

601 FORMAT (6(F11.5) ) 
C600 CONTINUE 
C------------ saving list of D, and Qext for the 5 wavelengths for 
PlOLlIng. ..-- 


Ge imum. ne.i) GOTO 702 


DO 700 J=2,K1 


WieiEca oO ALLO(I,2, 1), ALLO(J,4,1),ALLO(0,4,2) ,ALLG(0,; 4,3), ALLO(J 


1,4,4),ALLO(J,4,5) 


700 CONTINUE 
7O2 CONTINUE 
C>$2-+S Se Saving list of D32, and Qbar for the 5 wavelengths for 
DLOttinger 
cG SSS ES SF ES SS TS TS 8 TS 8 TS FE FE FE FE FE FE FF FF 5 FF FF 5 FF FF FF FF FE FF FS 5 FF FF FF 5 FE FS FF FF FF SF FS SE SE SE SS SE SS SS SS SS SS SS SS SS SE SES SS 
e: EVALUATING THE TRANSMISSION LOG RATIOS AND SSD's 
e 
‘e: SMALLEST=1000. 
K=ND+1 
DO 650 N=2,K 
RAT51(N)=ALLO(N, 4,5) /ALLQ(N, 4,1) 
RAT52 (N) =ALLQ(N, 4,5) /ALLQ(N, 4, 2) 
RAT53(N)=ALLO(N, 4,5) /ALLOQ(N, 4,3) 
RAT54 (N) =ALLO(N, 4,5) /ALLQ (N, 4, 4) 
RAT41(N)=ALLQ(N, 4, 4) /ALLQ(N, 4,1) 
RAT 42 (N) =ALLOQ(N, 4, 4) /ALLO(N, 4, 2) 
RAT43 (N) =ALLO(N, 4, 4) /ALLQ(N, 4,3) 
RAT 31 (N)=ALLO(N, 4, 3) /ALLOQ(N, 4,1) 
RAT32 (N) =ALLO(N, 4,3) /ALLQ(N, 4, 2) 
RAT21 (N)=ALLO(N, 4,2) /ALLQ(N, 4,1) 
RAT51A(N)=(RAT51 (N) -ALOG (T5) /ALOG(T1)) **2 
RAT52A (N) =(RAT52 (N) -ALOG (T5) /ALOG(T2) ) **2 
RAT53A (N) = (RAT53 (N) -ALOG (T5) /ALOG(T3) ) **2 
RAT5 4A (N) =(RAT54 (N) -ALOG(T5) /ALOG(T4) ) **2 
RAT41A(N)=(RAT41 (N) -ALOG(T4) /ALOG(T1) ) **2 
RAT 42A (N) =(RAT42 (N) -ALOG (T4) /ALOG (T2) ) **2 
RAT 4 3A (N) =(RAT43 (N) -ALOG (T4) /ALOG(T3) ) **2 
RAT 31A(N)=(RAT31 (N) -ALOG (T3) /ALOG(T1) ) **2 
RAT32A (N) =(RAT32 (N) -ALOG(T3) /ALOG(T2) ) **2 
RAT21A (N) =(RAT21 (N) -ALOG (T2) /ALOG(T1)) **2 
SSD(N) =RAT51A(N) +RAT52A(N) +RAT53A (N) +RAT54A(N) +RAT41A(N) +RAT4 2A (N) 
SSD(N) =SSD(N)+RAT43A(N) +RAT31A (N) +RAT32A (N) +RAT2Z1A(N) 
IF(SSD(N) .LT.SMALLEST) NBEST=N 
ITF (SSD(N).LT.SMALLEST) SMALLEST=SSD (N) 
650 CONTINUE 
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C sss========collecting the optimum results along all the best ones 
IF (SSD(NBEST) .LT.OPTSSD) OPTRDUST=RD@ST1 
IF (SSD(NBEST) .LT.OPTSSD) OPPS@6MA=SIGMA 
IF (SSD (NBEST).LT .@PESSD) OPTD32=A0EO(NBEST, 2.3 
TF (SSD (NBEST )2LT Seb tssb) OP fas P=sobD{( NSEs 1) 


c PRINTOUT OF TRANSMISSION RATIOS AND LOG QBAR RATIOS 
LE (NUMIONE* 1 )GO10 9 ie 
WRITE (2, 705) RDUST1, SIGMA, ALLO (NBEST, 2,1) ,SSD(NBEST) 


T51=ALOG (T5) /ALOG (T1) 
T52=ALOG (T5) /ALOG (T2) 
T53=ALOG (T5) /ALOG (T3) 
T54=ALOG (T5) /ALOG (T4) 
T41=ALOG (T4) /ALOG (T1) 
T42=ALOG (T4) /ALOG (T2) 
T43=ALOG (T4) /ALOG (T3) 
T31=ALOG (T3) /ALOG (T1) 
T32=ALOG (T3) /ALOG (T2) 
T21=ALOG (T2) /ALOG (T1) 
WROEE (2, /14)T51,RATSI(NBEST) , T51 
714 FORMAT (3(1x,F12.8)) 
WRITE (2, 7 4) 752, RATS2(NBEST)y eZ 
WREIE (2, (Payers, RATS3S (NEES) es 
WRITE (2,714)T54, RAT54 (NBEST) , T54 
WRITE (2, 714)T41,RAT41(NBEST),T41 
WRITE (2, 714)T42, RAT42 (NBEST) , T42 
WRITE (2, 7 i) 1.43;,RAT43 (NBEST War43 
WRITBE( 2,7) 14) 731, RAL3 IT (NBEST Mero 
WRITE (27714) 132,/RATS2Z(NBESL) i32 
WRITE (27714) T21,RATZI(NEEST) > 721 
v5 CONTINUE 


= ee recording all the best ones in file BestD32.dat ------ 
o Goren s3 3 

WRITE (>, 705) RDUST1, SIGMA, ALLO (NBBaw, 2, 1),SSD(NBEST) 
705 FORMAT (* m=",2f£6.4,'i Siogma="7i5na, 

1' gives best D32=',f6.4,' w/ SSD=',F12.7) 


e printing into screen SSD values to follow the run... 
WRITE (>; 707) SSD(NBEST) 


107 PORMAT (35%, "SSp— °, £12.7) 

383 continue 
RDUSTII=RDUST1II-SBBPIi 

6 CONTINUE 

C= 93334 ee eee ened. OOS nOne ROU ae —— 
RDUSTIR=RDUSTIR+STEPR 

8 CONTINUE 

SaaS ae —enee loop on RDUST I .=._—————__—-— 
SIGMA=SIGMA+STEPS 

9 CONTINUE 

C3 = 6 = oe = Se 2 8 SS = end ocp. on Sl cnaer=—————<_ === = ———— 

re 

C00. SSS6—+-—- SS recording and primgdng thesep@imum results -—---—— 


WRITE (*,801)WAV(1),T1,WAV(2),T2,WAV(3),T3,WAV(4),T4,WAV(5),T5 
WRETE (5, 801)WAV (1), TAWA), 12, WAV (3) oe, Wee), T4,WAV(5),T5S 
801 FORMAT (/,' for T(',£6.4,*)=",£6.4," T(";86e4, )—", tee 
Pe e's £6.4,")=", 26.4," TO aes, 9) =" FG od, (iG ae ee eee 
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WRITE (*, 802) OPTRDUST, OPTSIGMA, OPTD32, OPTSSD 
WRITE (5, 802) OPTRDUST, OPTSIGMA, OPTD32, OPTSSD 

610 2 Bonin tHe OPTIMUM RESULTS ARE:',/,'m=",216.4,. au 
Payvesedomomac—', i523," which gives D32=',16.4,' and =SSD= 3 
‘lite 


SUBROUTINE QAVG(M, NE, ND, SIGMA, NPTSZ,SD32, SDL32) 
BY KENNETH CASHDOLLAR, 1976 
QAVG CALCULATES D32 AND QBAR FOR A GIVEN QEXT VS D AND 
A GIVEN LOG-NORMAL SIZE DISTRIBUTION 


INPUT DATA = ND, SIGMA, D32, DEL32 


SOS) (ONO SA Ok Ske: 


ENEEGER M 

INTEGER ND 

REAL ALLQ, DQ, DIAM32 

COMMON ALLO (10000, 4,5) 
COMMON/ADQ/DQ (1000, 2) 
COMMON/AVG1/DIAM32 (1000) , WAV(5) 


DO Cie, DOC lL, Zy—OEXT 

NI=WAVELENGTH OR REFRACTIVE INDEX CHANGE 
ND IS NUMBER OF D32,QBAR CALCULATIONS 
SIGMA IS THE GEOMETRIC STANDARD DEVIATION 
D32 IS THE SURFACE WEIGHTED MEAN DIAMETER 
DEL32 ISmEiNCREMENT FOR D32 
DIST IS SURFACE WEIGHTED LOG-NORMAL SIZE DISTRIBUTION 

M IS THE MAXIMUM NUMBER OF POINTS IN Q-ARRAY FOR MONODISPERSE 


PRINT 240 
240 FORMAT (8X, 'D32 
i Pie oe Dax, ox, NMBR', 6X, ‘DIOX', 7X," DS52x." 
ioe, | OBAR xo 1 ty) 


WRITE (*, 1000) WAV (NT) 
1000 FORMAT(' ENTERING SUBROUTINE QAVG FOR LAMBDA=',F7.4,' MICRON’ ) 


C102. OF OOO) EO Oreo) (ae) ie 


D32=SD32 
DEL32=SDL32 


Cc NUMBER OF POINTS FOR PLOTS 2 AND 4 IS NPTS2 
NPTS2=ND+1 
DO 17 NN=1,ND 


DMIN=D32/SIGMA**3. 
DMAX=D32*SIGMA**2.5 
AS=ALOG (SIGMA) **2. 
DG=EXP (ALOG (D32) -2.5*AS) 
QOSUM=0.0 

DSUM=0 .0 
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OM Oe Ome) (OPO Ol 10) 16) 


XSUM=0.0 
DHOLD=0 . 0 
DNSUM=0. 0 
XNSUM=0.0 
NMBR=0 

MM=M-2 

DO 13 N=1,MM 


TE (DQ(N, 1). DRPeDMIN) Gesete is 
LE (DO(N, 1)3.Gl.. Dia Ge TOs 
IF (NMBR.GT.0) GO TO 242 


242 CONTINUE 

DD=DQO (N+1, 1) -DQ(N, 1) 
AR=(ALOG (DO (N#15 17 DG) = 427, 2 0 -AS) 
DIST=EXP (-AA) *DO(N+1,1) /SORT (6.2832*AS) 
OSUM=QSUM+ (DIST*DQO(N+1,2) +DHOLD*DQ(N,2))*DD/2.0 
DSUM=DSUM+ (DIST*DQO(N+1,1)+DHOLD*DQ(N,1))*DD/2.0 
XSUM=XSUM+ (DIST+DHOLD) *DD/2.0 
XNSUM=XNSUM+ (DIST/DQ (N+1, 1) **2+DHOLD/DQ(N,1)**2)/2.0*DD 
DNSUM=DNSUM+ (DIST/DO (N+1,1)+DHOLD/DO(N,1))/2.0*DD 
XHOLD=D0O (N+1,1) 
DHOLD=DIST 
NMBR=NMBR+1 

13 CONTINUE 

15 Continue 
OBAR=0SUM/XSUM 
D32X=DSUM/XSUM 


D10X=DNSUM/XNSUM 
ALLO (NN+1, 2,NI) =D32X 
ALLO (NN+1, 4, NI) =QBAR 


DIAM32 (NN+1) =D32 


OBARCHK (NN+1,1,NI) =D32X 
OBARCHK (NN+1, 2, NI) =OBAR 


D32=D32+DEL32 
17 COnitanveE 
WRITE (*, 1001) WAV (NT) 
C1001 FORMAT(' Leaving Subroutine QOAVG for Lambda=',F7.4,' micron') 


RETURN 
END 


a mm mm mm ei as a ee i Se SS i ee 
— ee a si a i Te ai SS SS SS 


typical input file (mieinp.dat) should be like... 
0.8000 0.7000 0.6000 O=s@@0 70-4000 
0.3130 0.4600 0.6328 0.7200 0.3506 
04 .100 
1,0) ley. 
Pie 02 © 


MeOO-07050  Lo0000 
200 MOO ly UOC CUM Ge Ouse 
PZ Ome O10 OSs Cr as0 Ono 
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APPENDIX B. ROCKET MOTOR DRAWINGS 


Four basic pieces formed the rocket motor: injector head, chamber section, 
turbulence ring, and nozzle. In this appendix all four drawings are given in this same 
sequence. The injector head drawing was the same, for 8, 6, 5, 4, 3 and 2-fuel elements. 
This part was not water cooled. The injector head’s central orifice was meant for the 
oxygen flow and the peripherical orifices for the kerosene. Three water cooled chamber 
sections were used to provide the desired combustor length. The turbulence ring was a 
solid disk with internal diameter smaller than the chamber section internal diameter, and 
was mounted between the first and the second chamber sections. The nozzle was also 
water cooled and employed conical convergent and divergent sections. The rocket motor 
also utilized two extra rings, just like the turbulence ring except that the internal diameter 
was equal to the chamber section internal diameter. They had a radial, 1/16” diameter 
hole to provide mesterent of the chamber pressure. These two rings were mounted 
adjacent to the injector head and between the second and third chamber sections, one on 


each side of the turbulence ring. 
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INJECTOR HEAD DRAWING 
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2 ee jived 
Std. 8 eet: [eee [oe | 1 ae | AP | oF gS ee i ee ne ae ee ee ot). C1. vid 
| eee ee es ae, ee ee aie Soe Ce eRe ee F3 
1-19) ¢ Shee IVY a ° 3 | : ‘ j Pe a ai os SS Sh fence fy RCE Bk a eg IC Ss U Ha a an ins SE I all t aie aie et 
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Al eos id S . Fe ee a a 
| ne - ; ee : _ Poi _ é ; ae 
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tt a Ee ame 7 io io < : 
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APPENDIX C. PRESSURE VS. TIME TRACES 


The following pressure transducer and diode voltage traces are from the actual 
firings whose soot concentration results were plotted in Chapter V. All plots have a run 
number on the top. Variables plotted were: 

Po2Choke — the oxygen sonically choked throat upstream pressure 

Pfuel Venturi — the kerosene cavitating venturi upstream pressure 

Po2zLow — the oxygen line pressure immediately before entering the injector head 
Pfuellow — the kerosene line pressure immediately before entering the injector head 
Pchamber — two lines for chamber pressure (at the injector head and at mid-chamber) 


He-Ne Voltage — the diode voltage from the helium-neon laser beam 


Also given are the venturi orifice diameters used in the fuel (kerosene) line and the 
actual (O/F) ratio for each firing. It should be noticed that for some runs (those with too 
low pressure difference between PfuelVenturi and PfuelLow) the venturi calibration 
formula will not match the (O/F) shown. That was because in these cases the venturi did 
not actually cavitate and the calibration did not apply. For these cases specific extra fuel 
flow calibrations were done reproducing upstream and downstrem venturi pressures to 


get reliable mass flow rates. 
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APPENDIX D. WHITE LIGHT TRANSMISSIONS TRACES 


The following five-wavelength power peak white light traces represent the main 
(O/F) ratios whose soot concentration results were relatively repeatable, since most of 
these kind of traces are similar. 

Runs whose transmissions traces are shown are: 
5-41 -- O/F = 0.333, 5-39 -- O/F = 0.494, 4-01 -- O/F = 0.575, 


4-15 — O/F = 0.754, 3-02 -- O/F = 0.909, 3-05 -- O/F = 1.097. 
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